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Magnetization steps (MST's) from a strongly diluted antiferromagnet on the square lattice are 
discussed theoretically. Thermal equilibrium, at temperature T and magnetic field B, is assumed. 
Two specific cluster models with the largest and second-largest exchange constants, j' 1 ^ and J^ 2 ' , 
respectively, are considered in detail. In the J1-J2 model, J^ 1 ' is the nearest-neighbor (NN) exchange 
constant Ji, and j' 2 ' is the second-neighbor exchange constant J2. In the J1-J3 model, J^=Ji, 
and is the third-neighbor exchange constant J3. For these two models, all cluster types of 
sizes n c <5 are identified, and their statistics is expressed using perimeter polynomials for cluster 
types. All cluster types with sizes n c >l give rise to MST's. Calculated curves of the isothermal 
magnetization M as a function of B are given for widely different ratios J^/j' 1 '. Some of the 
information contained in these magnetization curves is conveyed more directly by the derivative 
curves, dM/dB versus B. The series of peaks in the derivative curve is called the "MST spectrum." 
This spectrum is much simpler when the ratio J^'/j' 1 ^ is very small. A detailed discussion of the 
exchange-bond structure and MST spectra for cluster models with a very small ratio / is 
given in the following paper. The "parent" models of the J1-J2 model are the Ji and J2 models. 
For the J1-J3 model the parent models are the Ji and J3 models. Relations between the J1-J2 and 
the J1-J3 models and their respective parent models are discussed. The three parent cluster models 
(Ji-, J2-, and J3-models) are identical, except for the neighbor associated with the one exchange 
constant J that is included in the model. Common properties of such "isomorphic" cluster models 
are discussed. 

PACS numbers: 05.50.+q, 75.50.Ee, 71.70.Gm, 75.10.Jm, 75.10.Nr, 75.60.Ej 



I. INTRODUCTION 

An earlier theoretical papeni (hereafter I) on magne- 
tization steps from a diluted Heisenberg antiferromagnet 
on the square lattice was largely devoted to the nearest- 
neighbor (NN) cluster model. To obtain magnetization 
steps (MST's), the NN exchange constant J\ must be an- 
tiferromagnetic (AF). The present paper discusses cluster 
models with two exchange constants: the largest, called 
and the second-largest, J^ 2 \ Both exchange con- 
stants are assumed to be AF. 

Some of the discussion in the present paper applies to 
any JW-JW cluster model, regardless of the neighbors 
associated with JW and J^ 2 \ However, the main focus is 
on the J1-J2 and J1-J3 models. In these models J^'=J\, 
and j( 2 > is either the second-neighbor exchange constant 
J2, or the third- neighbor exchange constant J3. 

As in I, the following assumptions are made: 1) Ther- 
mal equilibrium, at temperature T and magnetic field B, 
prevails. 2) All the magnetic ions are identical. 3) The 
cation sites form a square lattice, and the magnetic ions 
are randomly distributed over these sites. 4) Only a frac- 
tion x of all cations sites are occupied by magnetic ions. 
This fraction is well below the site percolation concen- 
tration x c for the relevant cluster model. 5) None of the 
magnetic interactions is anisotropic. Background mate- 



rial for the present paper may be found in I, in a recent 
review^ and in earlier papers^*^ 

This paper is organized as follows. The principal re- 
sults are presented in the main text. Supplementary ma- 
terial is relegated to appendices. The main text starts 
with a brief discussion of the crucial role of cluster "con- 
figurations" in the theory. All cluster types in the J\- 
J2 and J1-J3 models, subject to the restriction n c <5 on 
the cluster size n c , are then identified. The statistics for 
these cluster types is expressed using perimeter polyno- 
mials (PP's) for cluster types. As discussed in I, the PP's 
for cluster types are analogous to the conventional PP's 
for cluster sizes £ 

The contribution of any cluster type to the magneti- 
zation M is calculated by combining the results for the 
energy eigenvalues with those for the cluster statistics of 
that cluster type. The total magnetization M(T, B) is 
the sum of the contributions from all cluster types. In 
some respects the derivative curve, dM/dB versus B, is 
more informative than the magnetization curve, M ver- 
sus B. Examples of calculated magnetization and deriva- 
tive curves at constant T are given for widely different 
ratios JW/JM. 

The statistics for cluster types is independent of the 
spin S of the individual magnetic ions. However, the en- 
ergy eigenvalues, and therefore the magnetization curve, 
depend on S. Although much of the discussion is for any 
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5, all the numerical examples are for 5=5/2, which is 
the appropriate value for the Mn 2+ and Fe 3+ ions. Both 
theses ions are 5-state ions, and they usually have low 
crystalline anisotropy. Such ions arc useful for testing 
theories in which anisotropy is neglected. 

The structures of the magnetization and derivative 
curves are much simpler when the ratio / JW is "very 
small." Cluster models with such widely different mag- 
nitudes of j( 2 ) and JW are called "lopsided cluster mod- 
els." In addition to their interesting physics, lopsided 
models are useful because they apply to many materi- 
als. Lopsided cluster models are discussed in detail in 
the following paper £ 

II. CLUSTER CONFIGURATIONS 
A. Configurations 

The calculation of the magnetization MiT, B) in any 
cluster model requires: 1) the identification of all cluster 
types c in that model; 2) a calculation of the average mag- 
netic moment /j, c (T,B) per realization^ for each cluster 
type c; and 3) an evaluation of the probabilities of finding 
the various cluster types. 

Before carrying out the first and third of these tasks 
it is necessary to identify all the "cluster configurations" 
that exist in the particular cluster model. These cluster 
configurations are the fundamental building blocks of the 
theory, and also of the computer programs that are used 
to implement the theory. Cluster configurations were dis- 
cussed in Sec. IIIB of I. The discussion below brings out 
some new features. 



1. Cluster configurations 

A spin cluster consists of a finite number of exchange- 
coupled magnetic ions (spins) that occupy a set of cation 
sites. Spin clusters are considered to have the same con- 
figuration if and only if the sets of cation sites occupied 
by these clusters can be obtained from each other by 
symmetry operations of the space group of the cation 
structure. Each one such set of cation sites is called a 
"realization" of the configuration. 

The symmetry operations of the space group of the 
cation structure are the only symmetry operations con- 
sidered in the present work. They will be referred to, 
simply, as "symmetry operations." The symmetry op- 
erations that are relevant to the present work are the 
operations of the P4m space group of the square lattice, 
including the lattice translations. 

Realizations of the same configuration have the follow- 
ing important geometrical property. Starting from one 
realization, a rigid object can be constructed by joining 
all pairs of cation sites in that realization by straight-line 
segments. The straight-line segments may be viewed as 
"struts" that give the object its rigidity. Rigid objects 



constructed in this manner from different realizations of 
the same configuration, either have identical shapes or 
are chiral isomers (mirror images) of each other. Because 
of this geometrical property, configurations of clusters are 
sometimes viewed as the geometrical shapes of clusters. 



2. Cluster configurations of one specific cluster model 

A cluster model is specified by the set of exchange con- 
stants (the J's) that are included in the model. Any spin 
cluster in this model consists of spins that are coupled 
to each other, but not to other spins, by the J's of the 
model. Thus, any two spins in a cluster are connected 
by at least one continuous path of exchange bonds asso- 
ciated with this set of J's. No continuous path of such 
exchange bonds is allowed to exist between spins in dif- 
ferent clusters. 

The restriction on the allowed J's is a restriction on the 
allowed sets of cation sites associates with the clusters of 
the model. For example, in the 3\-3i model, any cation 
site of a cluster must have a NN site or a 2nd-neighbor 
site in the same cluster. The cluster configurations of the 
model arc all the configurations of the sets of cation sites 
that are allowed by the J's of the model. 



3. Identical configurations in different cluster models 

Sometimes, several cluster models are considered. It 
may then happen that the set of J's in one cluster model 
is only a subset of all the exchange constants in another 
cluster model. Any cluster configuration in the cluster 
model with the fewer J's is then also a configuration in 
the model with the larger set of J's. For example, the 
exchange constant J\ of the NN cluster model is a subset 
of the exchange constants in the J\- J2 model. Any cluster 
configuration in the J\ model is also a configuration in 
the J1-J2 model. The converse is not true; many cluster 
configurations that exist in the J\-Ji model do not exist 
the Ji model. 

For the same reason, all cluster configurations that ex- 
ist in the J\ model also exist in the J1-J3 model. Fig- 
ure 2 of I shows an example of the same configuration in 
different cluster models. Figure 2(b) of I shows the con- 
figuration in the J\ model, and Figs. 2(c) and 2(d) show 
the same configuration in the J1-J2 and J1-J3 models, 
respectively. 



B. From cluster configurations to cluster types 

Once a cluster model is specified, it is necessary to 
identify the cluster types that exist in the model. In 
the present work this identification was carried out by a 
series of computer program. The first set of programs 
identified all the (different) configurations that exist in 
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the specified cluster model. For each configuration, a re- 
alization that has one spin at the origin was generated. 
This realization is called the "prototype" of the configu- 
ration. Only configurations with no more than 5 cation 
sites were considered explicitly in the present work. 

A cluster type c is specified by a cluster size n c and by 
a bond list A To identify the cluster types that exist in the 
model, the prototypes of all the (different) configurations 
were first classified by size, i.e., by the number of spins 
in the prototype. The next step was to generate the 
bond lists for all prototypes of a given size, n c . The final 
step was to identify all prototypes of the same size n c 
that have identical bond lists. Such prototypes are, by 
definition^ realizations of the same cluster type, c. In 
fact, they are the prototypes of all the configurations r c 
of cluster type c. 

To summarize, the classification of the prototypes of 
different configurations by both size and bond list leads 
to: 1) all cluster types c, for each cluster size n c ; 2) 
the bond list for each of these cluster types; and 3) the 
prototypes of all the configurations of each cluster type. 



In contrast to n rc , the perimeter v Ta depends not only 
on the configuration but also on the cluster model. If 
the same configuration exists in two cluster models, the 
perimeters in the two models are, in general, not equal. 
For example, any configuration that is present in the J\ 
model is also present in the J\-Ji model. The perimeter 
in the J\ model will be called the J\ perimeter. Given 
any realization of the configuration r c , the J\ perimeter 
is the number of (cation) sites that are NN's of the sites 
in the realization, but are not themselves sites of the 
realization. The perimeter in the J1-J2 model (called 
the J1-J2 perimeter) is the number of sites that are cither 
NN's and/or 2nd-neighbors of the sites in the realization, 
but are not themselves sites of the realization. Therefore, 
for the same configuration, the 3\-Ji perimeter is larger 
than the J\ perimeter. As a consequence, the probability 
that a randomly-chosen spin is in a realization of this 
configuration will be lower in the J1-J2 model than in 
the Ji model [See Eq. (4) of I] . Similar results apply to 
a configuration that exists in both the J\ and the J1-J3 
models. 



C. Statistics of cluster types 

The goal of the statistics is to find, for each cluster 
type c, the probability P c that a randomly-chosen spin 
is in one of the realizations of this cluster type. The 
main assumption is that the magnetic ions are randomly 
distributed over the cation sites. The procedure for cal- 
culating P c as a function of x was outlined in Sec. Ill C 
of I. The procedure starts from the configurations r c of 
cluster type c. 

Any realization of cluster type c must also be a re- 
alization of one of the configurations r c of that cluster 
type. The probability P rc that a randomly-chosen spin 
is in some realization of the configuration r c is given by 
Eq. (4) of I. This equation contains two parameters that 
depend on the configuration: the lattice-combinatorial 
parameter n rc , and the perimeter v rc . After these two 
parameters are evaluated for each of the configurations 
r c , the probability P c is obtained by summing P Ta over 
all the configurations r c of cluster type c. This sum is 
given by Eq. (5) of I. 

The lattice-combinatorial parameter n rc is the num- 
ber of (distinct) realizations of the configuration r c that 
have one spin at the origin. This n rc depends only on the 
configuration. If the same configuration exists in more 
than one cluster model, then the corresponding lattice- 
combinatorial parameters are the same in all these mod- 
els. The computer program that was used to obtain n Ta 
was based on the principle that all realizations of a con- 
figuration can be generated from the prototype of the 
configuration by applying the symmetry operations of 
the PAm space group, including the lattice translations. 
If the same realization was generated by different symme- 
try operations, the count for n rc included this realization 
only once. 



III. CLUSTER TYPES 

A. Generic and specific models 

A JW-JW cluster model in which the symmetry 
classes of the neighbors associated with and 
are not specified will be called a "generic" model. If 
the symmetry classes of the neighbors are specified, the 
cluster model is "specific" rather than generic. In this 
paper, only two specific jW-j( 2 ) models are considered: 
the J1-J2 model (with |Ji|>|J 2 |), and the J1-J3 model 
(with |Ji|>|J3|). The reason for focusing on these spe- 
cific models is that exchange constants often tend to de- 
crease as the distance r of the relevant neighbor increases. 
Although the decrease is not always monotonic^*^ we ex- 
pect that in the vast majority of diluted antiferromagnets 
with a square cation lattice, J^ 1 >=J\ and is either 
J2 or J 3 . 

The cluster types in the J1-J2 model are not trivially 
related to the cluster types in the J1-J3 model. That 
is, the bond lists for the two models cannot be obtained 
from each other by replacing the J2 bonds by J3 bonds. 
The site percolation concentrations for these two specific 
models are also different £ The non-trivial dependence of 
the bond lists, and hence of the cluster types, on the 
specific cluster model implies that bond lists and cluster 
types cannot be given for a generic model. They can only 
be given for a specific model. 

B. Parent cluster models 

The "parent" cluster models of the J1-J2 model are 
the Ji model and the J2 model, each of which has only 
one of exchange constants of the J1-J2 model. The J\- J2 
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FIG. 1: Cluster types of the J1-J2 model, up to cluster size n c =4. Solid circles represent spins. Solid and dotted lines represent 
Ji bonds and J2 bonds, respectively. The labels for the cluster types are discussed in the text. 
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FIG. 2: Cluster types of quintets (n c =5) in the J1-J2 model. 
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model is not a simple combination of its parent models. 
Similarly, the J1-J3 model is not a simple combination of 
the Ji and J3 models, which are its parent models. 

There are many interesting relations between the J\- 
Ji and J1-J3 models and their respective parent mod- 
els. To avoid repeated interruptions in the main text of 
the paper, discussions of these relations are relegated to 
Appendices. For the limited purpose of calculating mag- 
netization curves numerically, the results in these Ap- 
pendices arc not essential. However, these results give 
a deeper insight into the physics. Some of these results 
will be quoted, and used, in the main texts of the present 
paper and of the following paper. Appendix lAl describes 
the strong similarity (called "isomorphism" ) between the 
three parent models, i.e., the Ji, J2, and J3 models. 
These are the only "parent models" considered in the 
present work. 



C. Cluster types in the J1-J2 model 

1. Cluster types 

In the J1-J2 model on the square lattice, there are 67 
cluster types of sizes n c <5. These cluster types are shown 
in Figs, n and [21 Cluster types of the same size, n c , ap- 
pear in the same column. The four columns in Fig. ^ 
show the cluster types with n c =l,2,3, and 4. The sin- 
gle column in Fig. [2] shows the cluster types with n c =5. 
Spins are represented by solid circles, J\ bonds by solid 
lines, and Ji bonds by dotted lines. Only one config- 
uration for each cluster type is shown. The bond lists 
for these cluster types, which specify all intra-cluster ex- 
change interactions^ are given in Appendix 151 

Each cluster type is labeled by two numbers separated 
by a hyphen. The first number is the cluster size n c . The 
second is a serial number (SN) within this cluster size. 
Thus, the format for any label is [n c -(SN)]. There is only 
one type of single (type 1-1), but two types of pairs: 2-1 
which is a Ji-pair, and 2-2 which is a J2-pair. There are 
four types of triplets (n c =3), 15 types of quartets (n c =4), 
and 45 types of quintets (n c =5). 



2. Four categories of cluster types 

The number of cluster types in Figs. ^ and [2] is rather 
large. It is therefore useful to classify them. The clas- 
sification scheme is not unique. In one scheme the 67 
cluster types in Figs. ^ an d0 are divided into four broad 
categories: 

1. "single," which has no exchange bonds; 

2. "pure Ji" cluster types, with only J\ exchange 
bonds; 

3. "pure J2" cluster types, with only J2 bonds; 

4. "mixed" cluster types, with both Ji and Ji bonds. 



The single (type 1-1) is at the left of the bottom row 
of Fig. n The only pure Ji cluster types are the other 5 
cluster types in the same bottom row together with the 
four cluster types in the bottom row of Fig. [21 The only 
pure J\ cluster types are: 2-1, 3-1, 4-1, and 5-1. All the 
remaining 53 cluster types are "mixed" types. 

The pure J\ cluster types, and the pure Ji cluster 
types, are related to cluster types in the (parent) J\ and 
Ji models, respectively. These relations are discussed in 
Appendix [O 



D. Cluster types in the J1-J3 model 

In the J 1- J3 model there are 82 cluster types of sizes 
n c <5. They are shown in Figs. and The format is 
the same as in Figs.n]and[2] except that the dotted lines 
now represent J3 bonds. The bond lists for the cluster 
types in Figs. [3 and 0] are given in Appendix iBl 

The labels for the cluster types of the J1-J3 model have 
the same format as those for the cluster types in the J\- 
Ji model. Therefore, many of the labels used in the 
two models models are identical. When the same label 
appears in both models, it often refers to different cluster 
types (i.e., different bond lists). Unless it is clear from 
the context, it is then necessary to specify the cluster 
model to which the label refers. 

Once again the cluster types in Figs. [3] and 0] can be 
divided into four categories: 

1. the "single" (type 1-1), at the left of the bottom 
row of Fig. [3J 

2. the 9 "pure J3" cluster types, consisting of the 
other 5 cluster types in the bottom row of Fig. [3| 
together with the 4 cluster types in the bottom row 
of Fig. H 

3. the 5 "pure Ji" cluster types: 2-1, 3-2, 4-3, 4-5, 
and 5-5; 

4. the 67 remaining cluster types, which are all 
"mixed" types. 



IV. CLUSTER STATISTICS AND PERIMETER 
POLYNOMIALS 

A. Results for small clusters 

The probabilities P c as a function of x were obtained 
using the procedure discussed in Sec. Ill CI Only cluster 
types with n c <5 were considered. Some of the results 
for the J\-Ji model are shown Fig. [5] They include the 
P c 's for the single (type 1-1), for the two types of pairs 
(2-1 and 2-2), and for the four types of triplets (3-1, 3- 
2, 3-3 and 3-4). These labels for the cluster types refer 
to Fig. ^ Also shown in Fig. [S[ are: the sum P4 of the 
probabilities for all 15 quartet types in Fig. ^ the sum 
P5 of the probabilities for all the quintet types in Fig. [21 
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FIG. 3: Cluster types of the J1-J3 model, up to cluster size n c —4. The format, including the format of the labels for the cluster 
types, is similar to that in Fig. except that the dotted lines represent J3 bonds. 



and the sum P>5 of the probabilities for all cluster types 
with sizes n c >5. 

Figure shows the corresponding probabilities for the 
J1-J3 model. In this case the cluster types refer to Figs. 03 
and In both Fig. [S] and Fig. H3 the highest value of 
x is below the relevant site percolation concentration, 
x c =0.407 for the J\-J% model, and £ c =0.337 for the J\- 
J3 modeli2*i£ For the J\ model, x c =0.593. 



the J1-J2 model and the probabilities for the same clus- 
ter types in the (parent) ,J\ and J2 models. The relations 
between the probabilities P c for pure- Ji and pure- J3 clus- 
ter types in the J1-J3 model and the probabilities for the 
same cluster types in the (parent) J\ and J3 models are 
also discussed. 



THE MAGNETIZATION CURVE 



B. Perimeter polynomials 

As discussed in I, the probability P c can be expressed 
succinctly in the form 

Pc = n c x^- l D c {q), (1) 

where D c (q) is a polynomial in q—1— x, defined as the 
perimeter polynomial (PP) for cluster type c. Ap- 
pendix [51 gives the PP's for cluster types with sizes n c <5 
in the J1-J2 model. These are the cluster types in Figs.Q] 
and Appendix iBl also gives the PP's for cluster types 
of the J1-J3 model whose sizes are n c <5. These are the 
cluster types in Figs. and 0] 

For the J\- J2 model, the following check on the results 
for D c (q) in Appendix 151 was performed. The conven- 
tional PP's, D(q), for cluster sizes (not types) were given 
for this model by Peters et ali& As expected, the poly- 
nomial D(q) for any cluster size s is equal to the sum 

D c (q) of the PP's in Appendix iBlover all cluster types 
c that have the size n c =s. 

Appcndix[D]discusses some relations between the prob- 
abilities P c for the pure- Ji and purc-J2 cluster types in 



A. Calculation procedure 

The procedure for calculating the magnetization 
M(T, B) was discussed previouslyi^ Briefly, the Hamil- 
tonian of one realization 8 of each cluster type c is diag- 
onalized, and the results are used to obtain the average 
magnetic moment [i c (T,B) per realization. For singles, 
fi c (T,B) follows the Brillouin function for spin S. For 
each of the other cluster types, fJ, c (T, B) exhibits a series 
of MST's as a function of B at very low T. 

The total magnetization M(T, B) is a statistically- 
weighted sum of fJ. c (T, B). This sum is given by Eq. (13) 
of I, namely, 

M(T, B) = V N cPc (T, B) = A total V — Mc (T, B), 

n c 

c c 

(2) 

where N c =P c N to ta.\/n c is the population of cluster type 
c, and Atotai is the total number of spins. The quantities 
M, N c and A^otai are all either per unit mass or per unit 
volume. 

The infinite sum in Eq. @ cannot be evaluated exactly 
because /J, C (T, B) and P c (or N c ) are known only for a fi- 
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FIG. 4: Cluster types of quintets (n c =5) in the J1-J3 model. 



nite number of cluster types c. Usually they are known 
only for cluster types whose sizes n c are no larger than 
some maximum size n max . The sum in Eq. ( [2J) is there- 
fore truncated after the finite sum over cluster types with 
n c <n max is evaluated exactly. The remainder (REM) 
from clusters of sizes n c >n max , is then approximated by 
the "remainder correction" R(T,B). In the present work, 
n max =5, so that the REM is from clusters larger than 
quintets. The remainder correction R(T, B) was obtained 
by the corrective quintets (CQUIN's) method. The ap- 
plication of this method to the NN cluster model was 
discussed in I. For a model with two exchange constants 
the CQUIN's method is considerably more involved. 

The number of spins in the REM is P>5-/Vtotai- The ac- 
curacy of the CQUIN's method is not an important issue 



if P>5<Cl. Because P>5 increases with x, the accuracy is 
not a significant issue if x is sufficiently small. All mag- 
netization curves in the present paper are for z<0.09, 
and are based on the J1-J2 model. Under these condi- 
tions, -P>5<3.6%, so that errors in M resulting from the 
CQUIN's method are not significant. The description of 
the CQUIN's method is postponed to the following pa- 
per, which also includes some examples for higher x. 



B. The two reduced magnetic fields 

The two exchange constants, JW and j( 2 \ lead to 
two energy scales for the Zeeman energy. In analogy to 
Eq. (10) of I, the primary reduced magnetic field b\ is 
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FIG. 5: Probabilities P c as a function of x for some cluster types of the J\-J% model. The solid curves are for individual cluster 
types, labeled as in Fig. Q Curve 1-1 is for singles, curves 2-1 and 2-2 are for the two types of pairs, and curves 3-1 up to 3-4 
are for the four types of triplets. The dashed curves are probabilities for some combinations of cluster types: Pi is the sum of 
the P c 's for all the 15 quartet types in Fig. Q P5 is the sum of the P c 's for all 45 quintet types in Fig. and P>s is the sum 
of the probabilities for all cluster types of sizes n c >5. 




FIG. 6: Probabilities P c as a function of x for some cluster types of the J1-J3 model. The labels for the various curves are 
similar to those in Fig. EJ except that the cluster types refer to Figs. [3] and 0] for the J1-J3 model. 
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defined as 

bi=git B B/\jV>\. 
The secondary reduced magnetic field 62 is 

b 2 = gm B/\jW\. 
Thus, at any given B, 



h/b 2 = jW/jW 



< 1. 



(3a) 
(3b) 

(4) 



The reduced magnetization m is defined as in I. That is, 



m = M/M , 



(5) 



where Mo is the true saturation value of M. The re- 
duced parameters b\, 62, and to, will be used in plots 
and discussions of the magnetization curves. 



C. MST's from the two different types of pairs 

At a low T the calculated magnetization curve, M 
versus B, includes a superposition of many series of 
MST's. Each series arises from some cluster type c of 
size 2<n c <n max . The magnitude (AM) C of a magneti- 
zation jump at each MST in the series is proportional to 
the population N c of the relevant cluster type. For low 
x the largest jumps (AM) C are for j' 1 ) pairs and J' 2 ' 
pairs. In both J1-J2 and J1-J3 models these pairs are 
cluster types 2-1 and 2-2, respectively (see Figs. and 

EJi. 

The MST's from JW pairs occur at the primary re- 
duced fields 



6i = 2,4,6,...,4S. 



(6a) 



The MST's from pairs occur when the secondary 
reduced field has the same values, i.e., at 



2,4,6, 



AS. 



(6b) 



Experimental values of the magnetic fields B at the 
MST's from pairs are often used to determine JW and 
j( 2 ) _ The temperature requirement for resolving the 
MST's from pairs is k B T<\J^\. This is a more 
stringent requirement than &bT<| J*- 1 -* | for resolving the 
MST's from jW pairs. 

Equations l|6a|l and (|6b|) use the reduced fields b\ and 
62. Using the magnetic field B, instead of the reduced 
fields, the ranges of B for the MST series from the two 
types of pairs may or may not overlap. The condition for 
avoiding overlap is 



< 1/25. 



HI) 



(7) 



For magnetic ions with 5=5/2 (e.g., Mn 2+ or Fe 3+ ) over- 
lap is avoided if / jW <0.2. 




1.00 



0.95 - 



0.90 - 



0.85 



FIG. 7: Magnetization curve at T=0, calculated from the 
J1-J2 model using the parameters 2=0.01, S=5/2, and 
J2/Ji=0.28. The ordinate is the reduced magnetization, 
m=M/Mo, where Mq is the true saturation magnetiza- 
tion. The abscissa is the primary reduced magnetic field 
bi=g^BB /\ Ji|, up to 7.5. MST's from Ji pairs (cluster type 
2-1) are indicated by long upward arrows. MST's from J2 
pairs (cluster type 2-2) are indicated by shorter downward 
arrows. The inset shows the full magnetization curve, up to 
complete saturation. 



D. Two examples of magnetization curves for very 

low X 



The main purpose of the following two examples is 
to illustrate the dependence of the MST pattern on the 
ratio J^ 2 -* /J^\ The examples are for £=0.01 in the J\-J% 
cluster model. To optimize the resolution of the spectra, 
the examples are for T=0. The value 5=5/2 is assumed. 

Figures0and0are for J2/Ji=0.28. Figure0shows the 
reduced magnetization to as a function of the primary 
reduced magnetic field b\. MST's from J\ pairs (cluster 
type 2-1) are indicated by long arrows, and those from J2 
pairs (cluster type 2-2) by shorter arrows. The main part 
of Fig. shows the magnetization curve up to 61 = 7.5. 
The full magnetization curve is shown in the inset. 

Figure shows the derivative dM/dB, in normalized 
unitsjii in fields up to b\ = 7, corresponding to the main 
part of Fig. The upper abscissa scale is for the sec- 
ondary reduced magnetic field 62- The derivative peaks 
at MST's from Ji-pairs (type 2-1) and from J2-pairs 
(type 2-2), are indicated by long and short arrows, re- 
spectively. Their reduced magnetic fields arc given by 
Eqs. Il(:i;il) and (|6b|) . Because the ratio J2/Ji=0.28 is 
higher than 0.2, the field ranges for the series of MST's 
from the two types of pairs overlap. 

Figure [5] shows the magnetization (top) and derivative 
(bottom) curves for £=0.01 when the ratio J^j Ji=0.028, 
i.e., smaller by a factor of 10 compared to the ratio in 
Figs. [7] and [SJ All other parameters are the same as 
for Figs. and El Because the ratio J2/J1 is now well 
below 0.2, the MST series from J±- and J2-pairs occur in 
field ranges that do not overlap. In the "gap" between 
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FIG. 8: The derivative of the magnetization 
curve in Fig. |7| The ordinate is the derivative 
dm/d6i=(|Ji|/g/XBM )(dM/dS). The lower abscissa 
scale is for the primary reduced field b\. The upper abscissa 
scale is for the secondary reduced magnetic field 62. The 
derivative peaks at MST's from cluster types 2-1 (Ji pairs) 
and 2-2 (J2 pairs) are indicated by long and short arrows, 
respectively. 



these two field ranges, the magnetization (upper curve) 
exhibits a plateau of apparent saturation, labeled as m s . 
The apparent saturation value, m s =0.961, agrees with 
the value of m s obtained from the NN cluster model for 
this x (Ref[l|). However, in contrast to the NN cluster 
model, the apparent saturation is reached only after the 
series of MST's from the J2 pairs (type 2-2) is completed. 
Clearly, the MST pattern for J 2 /Ji=0.028 (Fig. is 
much simpler than that for J^j Ji=0.28 (Figs. and |HJ. 

Figure [7| for J2/ 'Ji=0.28 shows a short magnetization 
plateau immediately after the initial magnetization rise. 
This plateau, which ends at the first MST from J2 pairs 
(not J\ pairs), is not the plateau of apparent saturation 
predicted by the NN cluster model. The value m=0.925 
at the first short plateau in Fig.|7|is well below the appar- 
ent saturation value m s =0.961 in the NN cluster model. 
The reason for the lower value is that the series of MST's 
from J2 pairs has not been completed before the start of 
this short plateau. 



VI. THE MST SPECTRUM 

A. Spectrum 

The derivative dM/dB function of B exhibits a 
peak at each MST. The pattern of the peaks in the deriva- 
tive curve will be called the "MST spectrum." Figure |H| 
and the lower curve in Fig. [!|] are examples of such MST 
spectra. Each peak in dM/dB is a "spectral line." Two 
or more spectral lines associated with MST's arising from 
different cluster types may overlap. When a spectral line 
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FIG. 9: Zero-temperature magnetization curve (upper curve) 
and its derivative (lower curve) for x=0.01. These curves are 
calculated using the J1-J2 model and the ratio J2/J1 =0.028. 
As in Fig. |H| the lower and upper abscissa scales are for 61 
and 62, respectively. The left ordinate scale is for m. The 
right ordinate scale is for dm/d&i. The arrows indicate the 
locations of the MST's from Ji pairs (type 2-1). The MST's 
from J2 pairs (type 2-2) are bunched up at low fields. The 
plateau of apparent saturation is labeled as m a . 



is due to only one MST from one cluster type c, the in- 
tegral of the spectral line with respect to B is equal to 
the magnetization jump (AM) C . 

There are several advantages of using the spectrum. A 
plot of the spectrum, dM/dB versus B, is very effective 
in conveying information visually. Another advantage is 
that the calculated spectrum is the exact spectrum from 
cluster types with n c <n max . 

As discussed earlier, the infinite sum in Eq. J5J) is split 
into a sum over clusters with n c <n max , and a remainder 
(REM) from larger clusters (n c >5 in the present work). 
The reason for the split is that the Hamiltonians of clus- 
ters with n c >n max have not been diagonalized. The finite 
sum is evaluated exactly, but the REM is approximated 
by R(T, B). The derivative of R(T, B) does not give the 
spectral lines from the clusters with n c >n max . These 
lines can be obtained only if cluster Hamiltonians for 
n c >n max arc diagonalized. 

A wealth of information, such as accurate exchange 
constants, can be obtained from a comparison of an ex- 
perimental spectrum with the calculated exact spectrum 
for n c <n max . All the spectra shown in this paper are 
the derivatives of the exact finite sum up to n max = 5, 
without the remainder (see Ref . fl2|) . 



B. Additional examples of MST patterns and 
spectra 

The examples in Sec. IVDl for J 2 M=0.28 and 0.028, 
assumed £=0.01. The following examples are for the 
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FIG. 10: The reduced magnetization m (top curve) and the 
MST spectrum (lower curve) at T=0 for 2=0.09, S=5/2. 
These curves are for J2/ Ji=0.28. Only the results in the 
range &i<7.5 are shown. The cluster types responsible for 
some of the spectral lines are indicated. 



FIG. 11: Zero-temperature magnetization curve (top curve) 
and spectrum (lower curve) for £=0.09, calculated from the 
J1-J2 model when J2/J1 =0.028. The range of the primary 
reduce field is limited to &i<7.5. The plateau of apparent 
saturation is labeled as m 3 . 



same ratios Jij J\, but for higher x. All other parameters 
(5=5/2, T=0), and the cluster model (J1-J2), are the 
same as in Sec. IV Dl An example of a zero-temperature 
spectrum calculated from the J1-J3 model will be given 
in the following paper i£ Calculated spectra at finite tem- 
peratures will be shown in Ref. ^3 m connection with 
analysis of experimental data. 



1. Magnetization curves and spectra for £=0.09 

Magnetization curves and spectra for .t=0.09 are 
shown in Figs. EH and ^] From Fig. [5] the probabili- 
ties P c for all four triplet types increase when x changes 
from 0.01 to 0.09. The largest increase is for triplet type 

3- 3, which is a J\ pair attached to a third spin by a J2 
bond. As a result, spectral lines from type 3-3 triplets 
are readily seen in Figs. HOIandllll 

Figure ^] shows the MST pattern (top curve) and 
the MST spectrum (lower curve) for £=0.09 when 
J2/ Jl=0.28. The range of the primary reduced field is 
b\ <7.5. The cluster types responsible for some prominent 
spectral lines are indicated. Clearly, for J2/J1 = 0.28 
the spectrum is quite complicated because of the overlap 
between different series of MST's from different cluster- 
types. 

Figure 1111 shows the results for the same x when 
J2/J1 =0.028. For this much lower J2/J1 ratio, the spec- 
trum is much simpler. All the discernable spectral lines 
occur in two separate field ranges. The top of the low- 
field range is slightly above 61 =0.42 where the series from 
type 3-4 triplets ends. The bottom of the high-field range 
is near b\ =0.935 where a small MST from quartet type 

4- 2 is barely discernable. The two field ranges are sepa- 
rated by a gap, i.e., by a field range in which there are 



no discernable spectral lines. The absence of discernable 
lines implies that m has reached a plateau. The value 
?7!.=0.715 at this plateau agrees with the apparent satu- 
ration value m s calculated from the NN cluster model. 

In the field range of Fig. ^2 the J\ pairs, which exist 
both in the J\ and the J1-J2 models, give rise to large 
MST's at ^=2, 4, 6. These are the 2-1 lines in Fig. [TT] 
Near each 2-1 line there is also a line from the 3-3 triplets. 
The 3-3 lines do not exist in the J\ model. Each 3-3 line 
together with the stronger nearby 2-1 line may be viewed 
as a fine structure (FS) that has evolved from a single 
spectral line, due to J\ pairs, in the J\ model. The sep- 
aration A62, in the secondary reduced field, between the 
3-3 line and the nearby 2-1 line is of order 1. The corre- 
sponding magnetic field separation AB is g^/S.B^\ J 2 \. 

Figure HT1 also shows a spectral line at &i=7 labeled as 
3-2 & 3-1. It corresponds to the coincidence of the first 
MST from triplets of type 3-1 and the first MST from 
triplets of type 3-2. The "intensity" of this combined 
line is just the sum of the two intensities. The other 
lines from the 3-1 and 3-2 triplets, at 6i=9, 11, 13 and 
15 (above the field range of Fig. Ill|) . also coincide when 
72/^=0.028. FigurclU) shows that the 3-1 and 3-2 lines 
still coincide when Jij Ji=0.28. It can be shown that this 
remains true as long as J2/ </i<T. For additional details 
see Ref. Q 



2. Spectrum for x=0.20 

For £=0.20, approximately 40% of the spins are in clus- 
ters of sizes n c >5 (see Fig.|SJl. When such a large fraction 
of the spins are in the REM of Eq. J2J), the accuracy of the 
CQUIN's method of treating the magnetization from the 
REM is open to question. Under these circumstances, 
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FIG. 12: Zero-temperature spectrum in the J\-Jn model for 
J2/J1 =0.028 when £=0.20. Cluster types responsible for 
some spectral lines are indicated. 
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FIG. 13: Expanded view of the spectrum in Fig. E21 0=0.20) 
for the range 0<6i<2.5. 



a plot of the exact spectrum for clusters with n c <5 is 
probably preferable to a plot of the magnetization curve. 

The calculated spectrum for £=0.20 will be shown only 
for the low ratio J2/ Ji=0.028. FigurelT2*lshows this spec- 
trum in the range 0<&i<7.5, at T=0. Cluster types re- 
sponsible for some of the spectral lines arc indicated. An 
expanded view of the spectrum in the range &i<2.5 is 
shown in Fig. U3 The comparison of Figs. 1121 and 1131 
with Fig. ^3 shows that the increase of x from 0.09 to 
0.20 has led to the following changes: 

1. There are new disccrnablc lines. Some of these new 
lines are in the FS near &i=2. Two of the lines in 
this FS are from cluster type 4-12. One of these 
4-12 lines is slightly above the 2-1 line (from pure 
Ji pairs in the J1-J2 model), and and the other 4- 
12 line is slightly below it . Actually there is even 
a stronger new line in the same FS from quartets 
of type 4-10. However, this line is not resolved be- 



cause it is very close to the still stronger line from 
the 3-3 triplets. 
2. A gap, in which there are no discernable spectral 
lines, separates the low-field and high-field parts of 
the spectrum. On the scale of Fig. ED there are no 
discernable lines between &i=0.63 and &i=0.92. 



VII. LOPSIDED MODELS 

The spectra for a "very small" ratio J2/J1, are shown 
in Figs. El and lllllTSI These spectra are relatively simple. 
They suggest the following five features for such small 
J2/J1 ratios. 

1. The spectrum consists of a low- field part and a 
high- field part, separated by a gap in which there 
are no discernable spectral lines. 

2. In the field range of the gap, the magnetization ex- 
hibits apparent saturation, with an apparent sat- 
uration value m s equal to that given by the NN 
cluster model. 

3. In the high-field part of the spectrum, many spec- 
tral lines that exist in the parent J\ model develop 
a FS. The most conspicuous FS evolves from those 
lines in the J\ model that are due to Ji-pairs, i.e., 
the lines at &i=2, 4, 6, etc. 

4. Separations between adjacent lines in the FS that 
has evolved from a single line in the J\ model are of 
order A&2^1- The corresponding separations AB 
are of order g^BAB^ J 2 \. 

5. In the low- field part of the spectrum, the sepa- 
rations between adjacent lines are also of order 
A6 2 - 1, ot gf i B BAB~\J 2 \. 

The same five features are also found in simulations that 
use the J1-J3 model when the ratio J3/J1 is "very small." 

When the ratio J' 2 ' / is sufficiently small that 
the five features listed above appear in the spectrum, 
the jW-jt 1 ' model is called "lopsided." Spectra from 
lopsided models are discussed in detail in the follow- 
ing paper. Far from a mere curiosity, lopsided models 
actually apply to many materials. In Ref. Il'il which 
appears in this issue, the theoretical results for lop- 
sided models will be used to interpret data obtained in 
(C 3 NH 3 ) 2 Mn x Cdi_ x Cl4 near 20 mK. 
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APPENDIX A: ISOMORPHISM OF THE 
PARENT Ji, J 2 , AND J 3 MODELS 

Starting from Fig. 1 of I, and using symmetry 
arguments^ one can show that the cluster types of the 
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Ji-model, of the J2-model, and of the J3-model, are iden- 
tical except for the difference in the symmetry class of the 
neighbor associated with the (only one) exchange con- 
stant J that is included in the model. That is, the bond 
lists for all cluster types in the J2 model can be obtained 
from those in the J\ model by replacing the label 1 for 
NN's by the label 2 for 2nd neighbors. The bond lists 
in the J3 model are also the same, except that the label 
3 for the 3rd neighbors replaces the label 1. Different 
cluster models, each with only one exchange constant, 
whose bond lists are related in this manner will be called 
"isomorphic." Corresponding cluster types in isomorphic 
cluster models will be called isomorphic cluster types. 

The exchange part of the cluster Hamiltonian is speci- 
fied by the bond list for that cluster typei Cluster Hamil- 
tonians of isomorphic cluster types are identical except 
for the numerical value of the only J (see Ref . Il7l) . Each 
cluster type, except the single, leads to a series of MST's. 
The magnetic fields at the MST's originating from iso- 
morphic cluster types are the same except for a scale 
factor which is proportional to the J in the model. 

For the three parent cluster models considered here, 
the probabilities P c for isomorphic cluster types are the 
same. The site percolation concentrations for these three 
models are also the sained It is noteworthy that for the 
square lattice the J4-model is not isomorphic to the J±, 
J2, and J3 models. This difference can be understood 
from Fig. 1 of I. The number of 4th neighbors surrounding 
the central cation site is 8, compared to 4 for the 1st, 2nd, 
and 3rd neighbors. 



APPENDIX B: BOND LISTS AND PERIMETER 
POLYNOMIALS FOR THE Ji-J 2 AND THE J1-J3 
CLUSTER MODELS 

Tablcs[I]and[n]in give the bond lists, and the perimeter 
Polynomials, D c (q), for the cluster types of the J1-J2 
model. The cluster types c, limited to sizes n c <5, are 
those shown in Figs. ^ and [21 

Tables ITTT1 and llVl in this Appendix give the bond lists, 
and the perimeter Polynomials, _D c (g), for the cluster 
types of the J1-J3 model. The cluster types c are those 
shown in Figs. [3] and 



APPENDIX C: RELATIONS BETWEEN "PURE" 
CLUSTER TYPES AND CLUSTER TYPES IN 
THE PARENT MODELS 

Consider first the cluster types of the J1-J2 model 
(Figs. H an d |2J. At the bottom of these figures are the 
"single" (type 1-1) and the nine pure J2 cluster types. 
Together, they are identical to the 10 cluster types of 
the (parent) J2 model. As pointed out in Appendix ^ 
the cluster types of the J2 model are isomorphic to the 
cluster types of the J\ model. The pictorial representa- 
tions of the cluster types of the J\ model, shown earlier 
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TABLE I: Bond lists and perimeter polynomials D c (q) for 
cluster types of sizes l<n c <4 in the J\-J% model. The labels 
for the cluster types are as in Fig. Q 



in Fig. 3 of I, are therefore also applicable to the cluster 
types of the J2 model. 

There are only four pure J\ cluster types in the J\- 
J2 model. Together with the single, they comprise only 
a subset of the 10 cluster types of the parent J\ model 
(Fig. 3 of I). The reason why some cluster types of the 
J\ model are not among the pure J\ cluster types of the 
J\- J-2 model is the following. As discussed in Sec. Ill A 31 
all cluster configurations that are present in the J\ model 
are also present in the J1-J2 model. Some of these con- 
figurations have two J\ bonds that are: 1) connected to 
the same spin, and 2) make a 90° angle with each other. 
In the J 1- J2 model, each such configuration has at least 
one J2 bond, so that it is a configuration of a mixed clus- 
ter of the J1-J2 model, not a configuration of a pure Ji 
cluster. If all configurations of a particular cluster type c 
of the Ji model contain consecutive J\ bonds that make 
a 90° angle, then this cluster type is not among the pure 
Ji cluster types of the J1-J2 model. 

Two 3\ bonds that are connected to the same spin are 
either at a 90° angle or at a 180° angle with each other. 
Consider the pure J\ cluster types that exist in the J\-3% 
model. Each such cluster type has only one configuration 
in this model. All the spins in this configuration are on a 
single straight line. In the parent J\ model, on the other 
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{1110;201;20;0} 


8q 16 


5-29 


{1120;200;00;2} 


8q 19 +4q 20 


5-8 


{1110;221;00;0} 


4g 16 


5-30 


{1022;200;00;0} 


8q 21 


5-9 


{1112;201;20;0} 


8q 16 


5-31 


{1200;002;20;0} 


8g 20 +24q 21 + 16g 22 


5-10 


{1121;212;10;0} 


8q 14 


5-32 


{1200;000;20;2} 


16q 20 +32g 21 +8g 22 


5-11 


{1111;022;22;0} 


q 1G 


5-33 


{1200;000;22;0} 


2Aq 2Q 


5-12 


{1100;010;02;0} 


8q 18 


5-34 


{1220;000;02;0} 


8q 20 +8q 21 


5-13 


{1102;210;00;0} 


8q 18 


5-35 


{1022;000;22;0} 


4q 19 


5-14 


{1120;010;02;0} 


16g 18 


5-36 


{1010;200;01;0} 


16g 18 


5-15 


{1210;100;02;0} 


8q 17 +8q 1S 


5-37 


{1100;220;00;1} 


16g 17 +16g 18 


5-16 


{1120;210;02;0} 


8q 17 +8q ls 


5-38 


{2110;001;20;0} 


8g 18 


5-17 


{1122;210;00;0} 


8q 17 


5-39 


{0210;102;00;0} 


32g 20 


5-18 


{1110;220;02;0} 


8q 17 +8q 18 


5-40 


{2102;010;00;0} 


16g 19 


5-19 


{1110;222;00;0} 


8q 17 


5-41 


{0102;012;00;0} 


4g 18 + 12g 19 +8g 20 


5-20 


{2112;110;20;0} 


4g 16 


5-42 


{0220;202;00;0} 


4g 21 + 16g 22 + 12g 23 +2g 24 


5-21 


{1100;020;02;0} 


8q 20 


5-43 


{2022;200;00;0} 


8g 21 +12g 22 


5-22 


{1100;022;00;0} 


4g 19 


5-44 


{0222;220;00;0} 


8g 20 








5-45 


{2222;000;00;0} 


I 20 



TABLE II: Bond lists and perimeter polynomials, D c (q), for the quintet types (n c =5) of the J1-J2 model. The labels for the 
cluster types are as in Fig. 



hand, the same pure J\ cluster type, with the exception 
of cluster type 2-1, has several configurations. For exam- 
ple, cluster type 3 in Fig. 3 of I has two configurations in 
the Ji model (configurations 3a and 3/3 in Fig. 4 of I). 
The same cluster type exists in the J1-J2 model as clus- 
ter type 3-1 (Fig. of the present paper), but the only 
configuration is 3a, with the two J\ bonds at 180° angle. 

Similar results are obtained for the pure cluster types 
in the J1-J3 model: 1) The pure J3 cluster types of this 
model, together with the single, are identical to the 10 
cluster types of the parent J3 model. 2) The pure J\ 
cluster types, together with the single, are only a sub- 
set of the cluster types of the parent Ji model. In this 
case, the reason why they are only a subset is that in 
the J1-J3 model any configuration that has two consec- 
utive Ji bonds at 180° angle, also has at least one J3 
bond. Therefore, in the J1-J3 model it is a configuration 
of a mixed cluster, not of a pure J\ cluster. If all con- 
figurations of a particular cluster type c of the J\ model 
contain consecutive J\ bonds that make 180° angle, then 
this cluster type is not among the pure J\ cluster types 
of the J\- J3 model. 

To follow up on the previous example of cluster type 
3 in Fig. 3 of I, this cluster type is identical to cluster 
type 3-2 in the J1-J3 model (Fig.[3]of the present paper), 



which is a pure-Ji cluster type. Cluster type 3-2 of the 
J1-J3 model has only one configuration, namely, the 3/3 
configuration in which the two consecutive J± bonds are 
at 90° angle. 



APPENDIX D: COMMENTS ON THE 
PROBABILITIES OF PURE CLUSTER TYPES IN 
THE J1-J2 AND J1-J3 MODELS 

As discussed in Appendix [CI the pure J2 cluster types 
of the J1-J2 model are identical to cluster types of the 
parent Ji model. The configurations for each of these 
cluster types are also the same in the two models. How- 
ever, the probability P c for each of these cluster types 
is lower in the J\-Ji model than in the J2 model. The 
reason is that the J1-J2 perimeter for any configuration 
is larger than the J2 perimeter for the same configuration 
(see Sec. Ill Cjl . Similar remarks also apply to the proba- 
bilities for the pure J3 cluster types in the J1-J3 model, 
compared to the probabilities for the same cluster types 
in the parent J3 model. 

The case of the pure J\ cluster types is somewhat dif- 
ferent. Each of the pure J\ cluster types in the J1-J2 
model is identical to one of the cluster types in the par- 
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Cluster type, c 


Bond List 




1-1 


n 

I J 




2-1 


Hi 


2o 12 


2-2 


13) 


2« 13 


3-1 


111-31 


2g 16 


3-2 


1110) 


4g 15 


3-3 


113-01 


8g 16 +4g 17 


3-4 


13301 


4g 17 +2g 18 


4-1 


{113;31;0} 


2g 20 


4-2 


1 131-10-01 




4-3 


{110;01;0} 




4-4 


{111;30;0} 


4g 18 


4-5 


loii'iroi 


a 16 


4-6 


{110;33;0} 


8g 19 +4g 21 


4-7 


{113;30;0} 


4g 19 


4-8 


lllO'03-Ol 


16c 19 +8o 20 


4-9 


1113-00-0) 


8o 19 


4-10 


{130;03;0} 


4a 20 +8a 21 +2a 22 


4-11 


{130;00;3} 


16g 20 +16g 21 -|-4g 22 


4-12 


{133;00;0} 


12q 20 


4-13 


11 30-03-1 1 


2g 18 


4-14 


{310;01;0} 


4g 19 +8g 20 +2g 21 


4-15 


{330;03;0} 


4g 21 +8g 22 +2g 23 


4-16 


{333;00;0} 


4g 21 


4-17 


{033;33;0} 


3 20 



cnt Ji model. However, the probability P c is lower in 
the J1-J2 model, for two reasons. One is that, with the 
exception of cluster type 2-1, the number of configura- 
tions in the 3\-Ji model is smaller than in the parent 
Ji model (see Appendix IU|) . The second reason is that 
for those configurations that exist in both models, the 
J1-J2 perimeter is larger than the J\ perimeter. A larger 
perimeter implies a lower probability for the configura- 
tion. For the same two reasons the probability P c for 
any pure Ji cluster types in the J1-J3 model is lower 
than that for the same cluster type in the J x model. 



TABLE III: Bond lists and perimeter polynomials, D c (g), for 
cluster types of the J1-J3 model that have sizes l<n c <4. The 
labels for the cluster types are as in Fig. [3 
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Cluster Type, c 


Bond List 


D c (q) 


c 


Bond List 


D c (q) 


5-1 


{1133;310;01;0} 


2g 24 


5-30 


{1130;303;00;0} 


8g 22 +8g 24 


5-2 


{1310;130;01;0} 


8g 22 


5-31 


{1133;300;00;0} 


2q 22 


5-3 


{1130;010;01;0} 


8g 21 


5-32 


{1100;030;00;3} 


8q 22 +24g 23 +32g 24 +8g 25 


5-4 


{1133;010;01;0} 


4q 20 


5-33 


{1100;030;03;0} 


4g 22 +20g 23 +8g 24 +4g 25 


5-5 


{0110;101;00;0} 


4q 20 


5-34 


{1100;033;00;0} 


24g 23 


5-6 


{3110;101;00;3} 


4g 20 


5-35 


{1130;000;00;3} 


16g 23 +8g 24 


5-7 


{3110;101;00;0} 


4g 20 


5-36 


{1130;003;00;0} 


24g 23 +16g 24 


5-8 


{1131;310;00;0} 


8q 21 


5-37 


{1133;000;00;0} 


4g 22 


5-9 


{1311;100;00;3} 


4q 20 


5-38 


{1300;000;30;3} 


16g 24 +40g 25 +24g 26 +4g 27 


5-10 


{1110;301;00;0} 


8q 20 


5-39 


{1300;000;33;0} 


24g 24 +12g 25 


5-11 


{1101;013;10;0} 


8q 19 


5-40 


{1300;003;30;0} 


32g 24 +24g 25 +24g 26 +4g 27 


5-12 


{1111;300;00;3} 


q 2 ° 


5-41 


{1033;300;00;0} 


8g 24 +12g 25 


5-13 


{1130;310;03;0} 


8g 23 +4g 25 


5-42 


{1330;000;03;0} 


32g 24 +24g 25 


5-14 


{1133;310;00;0} 


8q 22 


5-43 


{1333;000;00;0} 


4g 23 


5-15 


{1310;100;03;0} 


16g 21 +8g 23 


5-44 


{1033;000;33;0} 


8g 23 


5-16 


{1130;010;03;0} 


8g 22 +8g 23 


5-45 


{1100;330;00;1} 


8g 22 +8g 23 +8g 24 +4g 25 


5-17 


{1103;310;00;0} 


16g 21 


5-46 


{3110;001;30;0} 


4g 23 


5-18 


{1313;100;00;0} 


8q 21 +8q 22 


5-47 


{1130;303;00;1} 


8g 21 


5-19 


{1100;010;03;0} 


8q 21 +8q 22 +8q 23 


5-48 


{1010;300;01;0} 


8g 21 +24g 22 +16g 23 +8g 24 


5-20 


{1103;010;00;0} 


16g 22 


5-49 


{3011;100;00;0} 


8g 22 +8g 23 


5-21 


{1110;303;00;0} 


16g 21 +8g 23 


5-50 


{1301;030;10;0} 


8g 21 


5-22 


{1110;003;30;0} 


8g 22 +4g 23 


5-51 


{0103;013;00;0} 


4g 22 +20g 23 +16g 24 + 12g 25 +2g 26 


5-23 


{1113;300;00;0} 


4g 21 


5-52 


{0310;103;00;0} 


8g 23 +32g 24 +24g 25 +4g 26 


5-24 


{0113;110;00;0} 


8q 20 


5-53 


{1303;030;10;0} 


16g 22 


5-25 


{1100;330;03;0} 


4g 22 +8g 24 +2g 26 


5-54 


{3103;010;00;0} 


32g 23 +24g 24 


5-26 


{1100;330;00;3} 


8g 23 +8g 24 +8g 25 +4g 26 


5-55 


{0330;303;00;0} 


4g 25 + 16g 26 +12g 27 +2g 28 


5-27 


{1100;333;00;0} 


4g 22 +8g 23 


5-56 


{3033;300;00;0} 


8g 25 +12g 26 


5-28 


{1100;330;03;3} 


4q 22 


5-57 


{0333;330;00;0} 


8g 24 


5-29 


{1130;300;00;3} 


8g 23 +4g 24 


5-58 


{3333;000;00;0} 


q 24 



TABLE IV: Bond lists and perimeter polynomials, D c (g), for quintet types (n c =5) of the J1-J3 model. The labels for the 
cluster types are as in Fig. 0] 



3-2 triplets produce additional lines at b\< 7. The values 
of 61 at these additonal lines depend on the ratio J2/ Ji . 
Conclusions based on numerical simulations for specific 
values of S and specific ratios J2/J1 fall short of being 
proven. However, the physical arguments in the following 
paper give strong additional support to these conclusions. 
All lattice sites of the square lattice in Fig. 1 of I can 
be reached from the central site by successive jumps be- 
tween NN sites. The sublattice consisting of all sites that 
can be reached from the central site by successive jumps 
between 2nd neighbors is, by itself, also a square lattice. 
Compared to the original square lattice, this "new" square 
lattice has a lattice constant that is a factor of v2 larger, 
and it is rotated by 45°. The J2 exchange bonds in the 
original square lattice are the NN exchange bonds in the 
new square lattice. All configurations, cluster types, clus- 



ter Hamiltonians, and cluster statistics for the original J2 
model are identical to those of the " J\ model" for the new 
square lattice. A similar argument applies to the sublattice 
consisting of all sites that can be reached from the central 
site by successive jumps between 3rd neighbors. 
In a given material the numerical value of J is specified by 
the symmetry class of the neighbor associated with J, i.e., 
it depends on whether J is J%, or J2, or J3. However, if 
the material is not specified then the numerical value of J 
is not known even if the symmetry class is known. In that 
case, J is just a parameter in the exchange Hamiltonian. As 
long as it remains just a parameter, there is no difference 
between the exchange Hamiltonians of isomorphic cluster 
types. 



